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ABSTRACT 

We present the stellar population and velocity dispersion gradients for a sample of 24 
brightest cluster galaxies (BCGs) in the nearby Universe for which we have obtained 
high quality long-slit spectra at the Gemini telescopes. With the aim of studying the 
possible connection between the formation of the BCGs and their host clusters, we ex- 
plore the relations between the stellar population gradients and properties of the host 
clusters as well as the possible connections between the stellar population gradients 
and other properties of the galaxies. We find mean stellar population gradients (nega- 
tive A[Z/H]/Zo<7 r gradient of -0.285±0.064; small positive Alog (age) /log r gradient 
of 0.069±0.049; and null A[E/Fe]/Zo# r gradient of -0.008±0.032) that are consistent 
with those of normal massive elliptical galaxies. However, we find a trend between 
metallicity gradients and velocity dispersion (with a negative slope of -1.616±0.539) 
that is not found for the most massive ellipticals. Furthermore, we find trends between 
the metallicity gradients and iGband luminosities (with a slope of 0.173±0.081) as 
well as the distance from the BCG to the X-ray peak of the host cluster (with a slope 
of -7.546±2.752). The latter indicates a possible relation between the formation of the 
cluster and that of the central galaxy. 

Key words: galaxies: evolution - galaxies: elliptical and lenticular, cD - galaxies: 
stellar content 



1 INTRODUCTION 

The formation of the supermassive galaxies at the centres of 
galaxy clusters, called brightest cluster galaxies (BCGs), is 
still one of the most challenging problems in galaxy forma- 
tion studies. The properties of these galaxies can be a result 
of their special location at the bottom of the cluster poten- 
tial well, or because BCGs occupy the massive end of the 
galaxy luminosity function (i.e. the properties can be driven 
by environment or mass). The comparison of the properties 
of BCGs with other ellipticals and with the properties of 
the cluster they reside in can help us elucidate which has 
the bigger influence on the evolution of these systems. 

While BCGs morphologically resemble elliptical galax- 
ies, their central surface brightnesses tend to be lower, and 
they follow a steeper Faber-Jackson relation (Faber & Jack- 
son 1976) than ordinary ellipticals (see, e.g. Von der Linden 
et al. 2007). They also have a luminosity function that dif- 
fers from the usual Schechter (1976) function that holds for 
normal cluster members (e.g. Hansen et al. 2005). The stars 
in BCGs have similar ages and metallicities than non-BCGs 
of the same mass, although some studies hint towards higher 
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[ct/Fe] measurements in BCGs (von der Linden et al. 2007; 
Loubser et al. 2009). 

Hierarchical models of galaxy formation predict that 
BCGs form through the (mainly late) assembling of small 
galaxies, and that the formation history of the BCG is 
closely linked to that of the host cluster (e.g., De Lucia & 
Blaizot 2007). Some observational results support this view: 
a significant alignment between the elongations of BCGs and 
their host clusters is observed in both the optical (Carter & 
Metcalfe 1980; Struble 1990; Plionis et al. 2003) and X-ray 
bands (Hashimoto, Henry & Boehringer 2008), and corre- 
lations between the BCG luminosity and cluster properties 
(for example X-ray temperature; Edge & Stewart 1991) have 
been reported. There is also a (weak) correlation between 
BCG mass and the mass of their host clusters which does 
not change significantly with redshift out to z ~ 0.8 (Edge 
1991; Collins & Mann 1998; Burke, Collins & Mann 2000; 
Brough et al. 2007; Stott et al. 2008; Whiley et al. 2008). 
Furthermore, an increase of BCG mass and size from z ~ 

0. 7 to z ~ 0.04 by a factor of 2 and 4, respectively, has 
been found (e.g. Bernardi 2009; Valentinuzzi et al. 2010), 
although compare with Whiley et al. (2008) and Collins et 
al. (2009) who did not find an increase in mass since z ~ 

1. However, it should be noted that the measurement of the 
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scale size of BCGs is notoriously difficult (Lauer et al. 2007; 
Stott et al. 2011 and references therein). 

Because the models predict that the majority of merg- 
ers happen at recent times, the gas content of the accreted 
galaxies is believed to be low (e.g. Dubinski 1998; Conroy, 
Wechsler & Kravtsov 2007; De Lucia & Blaizot 2007) and 
these mergers would therefore not change the central ages 
and metallicities of the BCGs. This would explain the lack 
of large differences between the stellar populations of BCGs 
and normal galaxies. However, these merger or accretion 
events are expected to change stellar population gradients, 
as dry minor mergers would deposit metal poor stars out- 
wards (Kawata et al. 2006). 

Stellar population gradients in ordinary early-type 
galaxies have been studied using line-strength indices by var- 
ious authors (Gorgas, Efstathiou & Aragon-Salamanca 1990; 
Gonzalez 1993; Fisher, Illingworth & Franx 1995a; Fisher, 
Franx & Illingworth 1995b; Cardiel, Gorgas & Aragon- 
Salamanca 1998a; Mehlert et al. 2003; Sanchez-Blazquez, 
Gorgas & Cardiel 2006b; Sanchez-Blazquez et al. 2007; and 
many more). A few studies have investigated the stellar pop- 
ulation gradients in small samples of BCGs (Gorgas et al. 
1990; Davidge & Grinder 1995 ; Fisher et al. 1995b; Carter, 
Bridges & Hau 1999; Mehlert et al. 2003; Sanchez-Blazquez 
et al. 2006b; Brough et al. 2007, Spolaor et al. 2009). Most of 
them either found no difference between the stellar popula- 
tion gradients of BCGs and those of the ordinary ellipticals 
(Mehlert et al. 2003), or slightly flatter metallicity gradients 
for the latter (e.g. Gorgas et al. 1990). 

However, a systematic study of the relations between 
the stellar population gradients and other properties of the 
cluster where they reside - e.g. the richness (parametrised 
with the X-ray luminosity) or the cluster mass - has not 
yet been carried out. This study presents a sample that is 
greater, in number, to all previous studies dealing with stel- 
lar population gradients in BCGs. Furthermore, our BCG 
sample has the advantage of having detailed X-ray values 
of the host cluster properties available from the literature, 
and significantly improved signal-to-noise (S/N) compared 
to previous samples (Crawford et al. 1999, and references 
therein). 

This paper is part of a series of papers investigating an 
overall sample of 49 BCGs in the nearby Universe for which 
we have obtained high S/N ratio, long-slit spectra on the 
Gemini and WHT telescopes. The spatially resolved kine- 
matics, central Single Stellar Population (SSP)-parameters, 
and the UV-upturn of the BCGs were presented in Loubser 
et al. 2008, Loubser et al. 2009, and Loubser & Sanchez- 
Blazquez 2011a (hereafter Paper 1, 2 and 3). The Mg2 gra- 
dients of a sub-sample of BCGs with sufficient S/N ratios 
were presented in Loubser & Sanchez-Blazquez 2011b (Pa- 
per 4). 

Here, we explore the possible correlation between the 
stellar population gradients and other properties of the clus- 
ters with the aim of clarifying if the galaxy gradients are 
shaped by the cluster assembly process. We fit SSP models, 
and present the stellar population gradients and their corre- 
lations with other properties, for a sub-sample of the BCGs 
with sufficient S/N ratios. We review the sample and data 
reduction, and discuss the data analysis in Section [2] The 
index and SSP gradients are derived and discussed in Sec- 
tions [3] and [4] respectively. We then correlate the findings 



with the host cluster properties in Section [6] and conclude 
our findings in Section [7| 



2 SAMPLE AND DATA REDUCTION 

The sample selection, spectroscopic data, and the reduction 
procedures were presented in Papers 1 and 2, and will not be 
repeated here. The procedure to derive the gradients from 
spatially binned spectra was presented in Paper 4, but will 
be briefly reviewed here for completeness. The galaxy spec- 
tra of the entire sample of 49 galaxies were binned in the 
spatial direction to a minimum S/N of 40 per A in the H/3 
region of the spectrum, ensuring a maximum error of ap- 
proximately 12 per cent on the measurement of this index. 
Twenty-six galaxies have four or more spatial bins when 
the central 0.5 arcsec, to each side, are excluded (thus in 
total the central 1.0 arcsec was excluded, comparable to 
the seeing). Those galaxies that did not meet this criterion 
were eliminated from the final sample. Thus, the gradients 
were investigated for 24 galaxies (2 were excluded because 
of emission contamination - see below). 

This subsample comprises of the dominant galaxies clos- 
est to the X-ray peaks in the centres of clusters and, for con- 
sistency, we call these galaxies BCGs to comply with recent 
literature (for example De Lucia & Blaizot 2007; Von der 
Linden et al. 2007; Brough et al. 2007Q 

Because some studies have found a possible difference 
between the evolution of the BCGs in high and low X-ray 
luminosity clusters (see Brough et al. 2002, 2005), it is use- 
ful to sample a range of cluster luminosities. Figure [I] shows 
the distribution of cluster velocity dispersion and the X-ray 
luminosities for a subsample (14 and 19 respectively) of the 
clusters studied here for which we have found data in the lit- 
erature. As can be seen, even in this reduced sample, the cov- 
erage in cluster velocity dispersion is quite large and range 
from o duster = 240 km s _1 (corresponding to galaxy groups) 
to o c i uat er = 1038 km s . We have 5 clusters with X-ray 
luminosity around or above the limit imposed by Brough 
ct al. (2002) to separate X-ray luminous and less-luminous 
clusters. These authors found that BCGs in X-ray luminous 
clusters have uniform absolute magnitudes (after correction 
for passive evolution) over redshifts 0.02 < z < 0.8, suggest- 
ing a lack of stellar mass evolution beyond that expected 
by passive evolution, while those BCGs in less X-ray lumi- 
nous clusters show significant scatter suggesting an increase 
in the mass up to a factor of ~ 4 since z ~ 1. According to 
this, most of our galaxies are in the range where a significant 
mass evolution is expected. 

The observations of this sub-sample were all obtained 
with Gemini using the multi-object spectrograph (GMOS) 
in the long-slit mode. Although we tried to place the slit 
along the major axis, in many cases (11 out of 24), there 
were no suitable bright guide stars in the area in front of 
the mask plane, which forced the slit to be rotated to an 
intermediate axis (see Table [T] for details). 

To use model predictions based on the Lick/IDS ab- 
sorption index system (Burnstein et al. 1984; Faber et al. 
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1985; Gorgas et al. 1993; Worthey et al. 1994), spectra 
were degraded to the wavelength dependent resolution of the 
Lick/IDS spectrograph where necessary (Worthey & Otta- 
viani 1997). Indices were also corrected for the broadening 
caused by the velocity dispersion of the galaxies (see the pro- 
cedures described in Paper 2). We measured line-strength 
indices from the flux-calibrated spectra and calculated the 
index errors according to the error equations presented in 
Cardiel et al. (1998b). We have also thoroughly investigated 
the possible influence of scattered light on the images by in- 
terpolating the counts measured across three unexposed re- 
gions in the GMOS 2D images. The scattered light signature 
is not necessarily spatially flat and was subtracted, with the 
normal background subtraction, from the images. We then 
compared the index measurements from the corrected and 
uncorrected images, and found that this contribution to the 
overall incident light is negligible within the spatial radii 
that we use here. 

We identified and corrected the spatially-resolved 
galaxy spectra contaminated by emission by using a com- 
bination of the ppxf (Cappellari & Emsellem 2004) and 
GANDALF (Sarzi et al. 2006) routine^ with the MILES stel- 
lar library (Sanchez-Blazquez et al. 2006c). We follow the 
same procedure as fully described in Paper 2. We found, 
and successfully removed, weak emission lines in NGC0541; 
NGC3311; NGC1713; NGC6173; and NGC7012. The emis- 
sion in ESO349-010 and MCG-02-12-039 could not be com- 
pletely removed without introducing erroneous features in 
some of the indices, and were excluded from further anal- 
ysis (a reliable Mg2 index gradient could be measured for 
MCG-02-12-039, but no H/3 gradient, and thus no stellar 
population parameters) . 

The half-light radii (r c ) were calculated from the 
2MASS catalogue, except ESO303-005 for which it was cal- 
culated by fitting a de Vaucouleurs profile to a cut in the 
spatial direction of a 2D spectrum, using the VAUCOUL task 
in the reduceme package]^] The half-light radii were com- 
puted from the 2MASS K-band 20th magnitude arcsec 2 
isophotal radius using the formula by Jarrett et al. (2003): 
log(r e ) ~ log(rK2o) ~ 0.4. For old stellar populations, these 
half-light radii do not differ much from those derived using 
the optical bands (Jarrett et al. 2003), and the typical errors 
on r e propagated from tk20 is less than 0.08 and therefore 
negligible]^] We projected this radius on to the major axis 
(MA), taking into account the ellipticity, as: 



1 - e | cos(\ PA - MA 



(1) 



Figure 1. The distribution of the host cluster velocity dispersions 
and X-ray luminosities (all measured in the 0.1 - 2.4 keV band; 
where known from the literature) . The vertical line in the second 
histogram represent the X-ray luminosity used in Brough et al. 
(2002) to separate X-ray luminous and less luminous clusters. For 
the references and other host cluster properties, see Table [5] 



2 We make use of the corresponding PPXF and GANDALF IDL 
(Interactive Data Language) codes which can be retrieved at 
http: /www . leidenuniv.nl/sauron/ 

° Keduceme is an astronomical data reduction package, specialis- 
ing in the analysis of long-slit spectroscopy data. It was developed 
by N. Cardiel and J. Gorgas (Cardiel 1999) ( htt p~/www ■ ucm ■ es/| 
inf o/Astrof /sof tware/reduceme/reduceme .html ). 
* However, sec Htott et al. (2011) for an extensive discussion on 
the uncertainties of measuring r e for BCGs. For our purpose, 
we only use the r e measurements to obtain an estimate of how 
far our measured gradients reach compared to the extend of the 
galaxy, and it therefore does not change the results of our gradient 
correlations. 
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with e the ellipticity (data from NED) , r c the radius contain- 
ing half the light of the galaxy (computed from the 2MASS 
if-band 20th magnitude arcsec -2 isophotal radius as de- 
scribed in Paper 1), PA the slit position axis, and MA the 
major axis. The fractions of a e which the radial profiles mea- 
sured in this work spans, are listed in Table [I] The central 
radial velocity and velocity dispersion values of the galax- 
ies were derived for regions with the size of o e /8 (see Paper 
1), and with the galaxy centres defined as the luminosity 
peaks. The index values at each radius were measured at the 
local velocity and velocity dispersion for that radius (aver- 
aged over the bin) as presented in Paper 1. The luminosity- 
weighted centres of the spatial bins were used to plot the 
gradients in this paper. 

Note that, despite our derived gradients reaching on 
average a fraction of 0.4 of the a e of the galaxy, Brough 
et al. (2007) showed that SSP gradients in BCGs do not 
deviate from a power law up to three times a e . Therefore, 
our results are valid even if very large radii are not reached. 



3 INDEX GRADIENTS 

Lick indices were measured according to the definitions of 
Trager et al. (1998). We plot the H/3, Mg b , Fe5270 and 
Fe5335 index measurements in magnitudes, obtained by us- 
ing 



I' = 2.5 log 



(AA C 



+ 1 



(2) 



where (AA C ) represents the width of the central bandpass, 
and I and I' are the index measurement in A and magni- 
tudes respectively. 

As a consistency test, we compare the index measure- 
ments of the central bins with the central measurements 
(in a c /8 apertures as presented in Paper 2) in Appendix [A| 
The index gradients of one of the galaxies (NGC4839) could 
also be compared to previous measurements in the litera- 
ture from Fisher et al. (1995b) and Mehlert et al. (2000), as 
shown in Figure |A2| The profiles compare well within the 
errors and small differences can be as a result of spectral 
resolution and slit width. 

We represent the index gradients as / = a + b\og(^ L ), 

where I is the index (in magnitudes) and a e is the effective 
radius as presented in Table[T] Gradients are taken to be sig- 
nificantly different from zero if the slope is greater than three 
times the la error on that gradient. Gradients (and the error 
on the gradient) were obtained by a linear least-squares fit 
using a Marquardt-Levenberg algorithm (Marquardt 1963) 
to the radial profiles. To take the errors on the individual 
values into account, statistical t-tests were also performed 
on all the slopes to assess if a real slope is present or if it 
is zero (as a null hypothesis). For our degrees of freedom, a 
t-value larger than 1.96 already gives a probability of lower 
than five percent that the correlation between the two vari- 
ables is by chance. P is the probability of being wrong in 
concluding that there is a true correlation (i.e. the probabil- 
ity of falsely rejecting the null hypothesis). 

We find only two galaxies with significant H/3 gradients 



MA was known) , we find a mean H/3 gradient and standard 
error of 0.001 ± 0.010. Similar results are found for normal 
elliptical galaxies (Mehlert et al. 2003; Sanchez-Blazquez et 
al. 2006b; Kuntschner et al. 2006, 2010). When we calculate 
the mean including all the derived H/3 gradients, regardless 
of slit orientation, we find 0.003 ± 0.012. The R/3 gradient 
of ESO488-027 showed an asymmetrical profile (see Section 

We find six galaxies with significant Mgb gradients, nine 
with significant Fe5270 gradients, and five with significant 
Fe5335 gradients (shown in Figures B3 to B6 in Appendix 



(shown in Figures Bl and B2 in Appendix [B] and Table 
[2j. When we separate the gradients of the galaxies where 
the slit was placed within 15 degrees of the MA (where the 



[B]and Table [2}. When we calculate the mean including all 
the derived Mgt, Fe5270 and Fe5335 gradients, regardless of 
slit orientation, we find -0.017 ± 0.022, -0.014 ± 0.015 and 
-0.010 ± 0.021 respectively. We again note asymmetrical 
profiles for some of the galaxies. 



4 SSP-EQUIVALENT PARAMETER 
GRADIENTS 

To calculate the ages, metallicities ([Z/H]), and a- 
enhancement ratios ([E/Fe]), we compare our derived line- 
strength indices with the predictions of Thomas, Maraston 
& Bender (2003) and Thomas, Maraston & Korn (2004) 
and using the Korn, Maraston & Thomas (2005) model at- 
mospheres, with a slightly modified method from the one 
presented by Trager et al. (2000). Variations of the indices 
with chemical partitions departing from solar were included, 
where the "E" group contains O, Ne, Mg, Si, S, Ar, Ca, 
Ti, Na and N. Ages, metallicities and a-enhancement ratios 
were derived for the BCG sample using the indices <Fc>Q 
H/3 and Mgb. Then, a \ 2 -minimisation was applied to find 
the combination of SSPs that best reproduced the four in- 
dices simultaneously. Errors on the parameters were calcu- 
lated by performing 50 Monte-Carlo simulations in which, 
each time, the indices were displaced by an amount given 
by a Gaussian probability distribution with a width equal 
to the errors on these indices. 

We take the stellar population gradients as the slope of 
a linear fit, inversely weighted by the errors, to the relation 
/ = a + 61og(r) where / is the SSP-equivalent parameter. 
This enables a direct comparison with previous studies, for 
example Brough et al. (2007). We fit the stellar population 
gradients using the same procedure as used for the index 
gradients, and once again exclude the central 0.5 arcsec (to 
each side) from the fit and take the gradient to be signif- 
icant if it is larger than three times the la error on that 
gradient. The SSP gradients are tabulated in Table [3] and 
shown in Figures |B7| to |B"T0| in Appendix [B| We only indi- 
cate the gradient fit in the plots (with a line) if it is bigger 
that 3 times the la error on the gradient. The age gradients 
of ESO146-028 and ESO303-005 showed asymmetrical pro- 
files. We have investigated these galaxies (with respect to 
their emission, slit placement, radial kinematic profiles (see 
paper 2), central stellar populations and host cluster X-ray 
properties) without finding any reason for the asymmetric 
gradients. Since these galaxies have few data points, the two 
galaxy sides were not fitted separately. 



<Fe> =(Fe5270 + Fe5335)/2, (Gonzalez 1993). 
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Table 1. Galaxies observed with the Gemini North and South telescopes (see Paper 1 for more detailed properties of the galaxies). The 
position angle (PA) is given as deg E of N. The galaxies marked with a * is where the slit PA was placed within 15 degrees from the 
major axis (MA; where the MA was known). Typical errors on r e is less than one per cent. The last column lists the fraction of the 
effective half-light radii spanned by the radial profiles measured in this work. ESO349-010 and MCG-02-12-039 are included in the table 
even though the emission lines could not be fully removed. 
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.77 


0.152 


0.024±0.012 


2.06 


0.109 


NGC0533 


0.035±0.014 


2.56 


0.043 


-0.014±0.013 


-1 


05 


0.334 


-0.001±0.004 


-0 


.12 


0.907 


0.007±0.014 


0.53 


0.614 


NGC0541 


-0.001±0.005 


-0.05 


0.960 


-0.004±0.007 


-0 


56 


0.599 


-0.001±0.003 


-0 


.20 


0.847 


0.008±0.005 


1.63 


0.164 


NGC1399 


-0-001±0.002 


-0.13 


0.901 


-0.045±0.002* 


-20 


61 


<0.0001 


-0.013±0.002* 


-8 


.50 


<0.0001 


-0.030±0.003* 


-11.26 


<0.0001 


NGC1713 


0.001±0.004 


0.28 


0.784 


-0.035±0.007* 


-5 


15 


0.002 


-0.018±0.005* 


-3 


.67 


0.003 


-0.019±0.004* 


-4.40 


0.001 


NGC2832 


0.001±0.006 


0.02 


0.983 


-0.065±0.017* 


-3 


76 


0.002 


-0.032±0.008* 


-4 


.24 


0.001 


-0.057±0.012* 


-4.75 


0.001 


NGC3311 


-0.015±0.004* 


-4.23 


0.008 


0.026±0.019 


1 


41 


0.217 


0.004±0.006 





.63 


0.556 


-0.006±0.005 


-1.07 


0.331 


NGC4839 


0.008±0.003 


3.09 


0.018 


-0.015±0.010 


-1 


47 


0.238 


-0.012±0.003* 


-3 


.74 


0.007 


-0.007±0.004 


-1.87 


0.104 


NGC6173 


-0.006±0.020 


-0.32 


0.780 


-0.014±0.027 


-0 


53 


0.651 


-0.014±0.006 


-2 


56 


0.124 


0.008±0.012 


0.64 


0.567 


NGC6269 


O.OlOiO.Oll 


0.98 


0.364 


-0.016±0.007 


-2 


22 


0.068 


-0.013±0.005 


-2 


59 


0.041 


-0.016±0.011 


-1.44 


0.201 


NGC7012 


0.007±0.007 


1.05 


0.309 


-0.025±0.007* 


-3 


59 


0.002 


-0.024±0.007* 


-3 


.39 


0.004 


-0.011±0.006 


-1.90 


0.076 


PGC004072 


0.011±0.009 


1.31 


0.237 


-0.022±0.005* 


-4 


14 


0.054 


-0.013±0.017 


-0 


.76 


0.528 


0.003±0.027 


0.11 


0.923 


PGC030223 


0.002±0.018 


0.08 


0.943 


-0.028±0.014 


-2 


07 


0.085 


-0.062±0.019* 


-3 


31 


0.016 


-0.041±0.007* 


-5.53 


0.002 


PGC072804 


-0.014±0.009 


-1.55 


0.181 


-0.048±0.019 


-2 




0.051 


0.001±0.002 





.59 


0.581 


0.006±0.007 


0.83 


0.442 


UGC02232 


-0.013±0.016 


-0.82 


0.456 


-0.013±0.016 


-0 


84 


0.449 


0.008±0.020 


0.38 


0.723 


-0.024±0.029 


-0.85 


0.444 


UGC05515 


-0.009±0.040 


-0.22 


0.845 


-0.024±0.020 


-1 


22 


0.346 


-0.030±0.009* 


-3.29 


0.081 


-0.013±0.004* 


-3.19 


0.086 



Mean 0.003 -0.017 -0.014 -0.010 

Std. dev 0.012 0.022 0.015 0.021 

Std. err 0.002 0.005 0.003 0.004 



Table 2. The BCG H/3, Mg;,, Fe5270 and Fc5335 gradients. A ★ indicates the cases where the gradient is bigger than three times the 
error on that gradient. 



Some galaxies (ESO552-020, IC1633, IC5358, 
NGC0533, NGC1399, NGC2832 - none of which show 
nebular emission) and NGC7012 (where the emission 
was removed) reached the upper age limit of the stellar 
population models used. Note that the oldest ages in certain 
models, such as the ones used here, are older than the cur- 



rent age of the Universe (see discussion in Sanchez-Blazquez 
et al. 2009). However, all the interpretations of this study 
are based on relative differences in ages which are much 
more reliable than absolute values. 
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Object log (ago) gradient t P [E/Fe] gradient t P [Z/H] gradient t P 



ESO146-028 


0.012±0.301 


0.04 


0.969 


0.261±0.155 


1.69 


0.126 


-0.107±0.166 


-0.65 


0.534 


ESO303-005 


0.587±0.598 


0.98 


0.381 


0.136±0.235 


0.58 


0.592 


-0.573±0.360 


-1.59 


0.187 


ESO488-027 


0.122±0.147 


0.83 


0.421 


-0.049±0.085 


-0.57 


0.577 


-0.226±0.081 


-2.77 


0.015 


ESO552-020 


-0.054±0.040 


-1.35 


0.202 


-0.003±0.049 


-0.07 


0.946 


-0.120±0.070 


-1.73 


0.110 


GSC555700266 


0.369±0.526 


0.70 


0.522 


0.044±0.146 


0.30 


0.781 


-0.267±0.418 


-0.64 


0.558 


IC1633 


0.094±0.043 


2.21 


0.031 


-0.061±0.031 


-1.98 


0.053 


-0.301±0.048* 


-6.23 


<0.0001 


IC4765 


-0.164±0.039* 


-4.21 


0.052 


-0.075±0.059 


-1.28 


0.329 


0.136±0.182 


0.75 


0.532 


IC5358 


-0.084±0.112 


-0.75 


0.467 


0.185±0.128 


1.45 


0.172 


-0.429±0.136* 


-3.17 


0.007 


Leda094683 


0.261±0.380 


0.69 


0.530 


-0.058±0.063 


-0.93 


0.405 


-0.761±0.493 


-1.54 


0.198 


NGC0533 


-0.471±0.196 


-2.41 


0.043 


-0.033±0.063 


-0.52 


0.615 


-0.141±0.071 


-1.98 


0.083 


NGC0541 


-0.014±0.128 


-0.11 


0.920 


-0.083±0.041 


-2.02 


0.100 


0.046±0.146 


0.32 


0.763 


NGC1399 


0.104±0.025* 


4.22 


<0.0001 


0.035±0.022 


1.62 


0.108 


-0.501±0.030* 


-16.51 


<0.0001 


NGC1713 


0.059±0.064 


0.93 


0.372 


0.003±0.053 


0.05 


0.959 


-0.475±0.080* 


-5.93 


<0.0001 


IN ZOOZ 


U.U44±U.U /V 


0.02 


0. 540 


n mnJ.n ncQ 

u. iuy±u.UDo 


2.08 


0.051 


— U.D /o±U.14y* 


—3.87 


0.01 


NGC3311 


0.255±0.128 


2.00 


0.102 


0.070±0.131 


0.54 


0.616 


-0.053±0.167 


-0.32 


0.764 


NGC4839 


-0.081±0.055 


-1.48 


0.181 


0.035±0.039 


0.89 


0.403 


-0.122±0.074 


-1.65 


0.144 


NGC6173 


0.094±0.138 


0.68 


0.532 


-0.196±0.067 


-2.91 


0.044 


-0.010±0.072 


-0.14 


0.898 


NGC6269 


-0.009±0.393 


-0.02 


0.983 


0.127±0.087 


1.45 


0.197 


-0.284±0.289 


-0.99 


0.363 


NGC7012 


0.020±0.112 


0.18 


0.858 


0.081±0.061 


1.33 


0.202 


-0.404±0.129* 


-3.13 


0.006 


PGC004072 


0.141±0.317 


0.44 


0.701 


-0.045±0.195 


-0.23 


0.839 


-0.306±0.165 


-1.85 


0.205 


PGC030223 


-0.166±0.273 


-0.61 


0.566 


-0.009±0.135 


-0.06 


0.951 


-0.725±0.185* 


-3.93 


0.008 


PGC072804 


0.638±0.208* 


3.06 


0.028 


-0.358±0.128 


-2.80 


0.038 


-0.685±0.215* 


-3.18 


0.025 


UGC02232 


-0.220±0.240 


-0.92 


0.428 


-0.442±0.036* 


-12.15 


0.001 


0.614±0.440 


1.40 


0.257 


UGC05515 


0.131±0.358 


0.36 


0.751 


0.123±0.069 


1.78 


0.216 


-0.578±0.001* 


-628.13 


<0.0001 


Mean 


0.069 






-0.008 






-0.285 






Std. dev 


0.241 






0.157 






0.315 






Std. err 


0.049 






0.032 






0.064 







Table 3. Age, mctallicity and ct-enhancement gradients for our sample of BCGs. A * indicates where the gradient is bigger than three 
times the error on that gradient. 



5 RESULTS 

In the following sections, we will investigate if the stellar 
population gradients of BCGs are shaped more by the galaxy 
properties or by the environment. To do this, we will explore 
the correlations between the stellar population gradients and 
the properties of both, the galaxies themselves, and their 
host clusters. 



5.1 Metallicity gradients and correlations with 
other galaxy properties 

The distribution of metallicity gradients for our sample of 
BCG galaxies is represented in Figure[5] while the numerical 
values together with the errors are listed in Table [3] It is 
immediately seen that our sample of galaxies show a large 
variety in metallicity gradients. This has already been shown 
in other studies of especially resolved stellar populations in 
BCGs (e.g. Gorgas et al. 1990; Brough et al. 2007), but 
we are showing that this relation also applies to a much 
larger sample here. The variety of the metallicity gradients 
for the BCG galaxies contrasts with their small dispersion 
in their aperture luminosities (Sandage 1972; Gunn & Oke 
1975; Hoessel & Schneider 1985; Postman & Lauer 1995). 

From Table [3] it can be seen that we find only nine 
significant metallicity gradients. When we include all the 
derived metallicity gradients, we find a mean gradient of 
-0.285±0.064 (with a standard deviation of 0.315), which 
compares very well to -0.31 ± 0.05 found by Brough et 
al. (2007). For normal ellipticals, Sanchez-Blazquez et al. 
(2007), Gorgas et al. (1990) and Fisher et al. (1995b) 
found mean [Z/H] gradients of A [Z/H] /log r = -0.306 ± 
0.133, -0.23 ±0.09 and -0.25 ±0.10, respectively, , that is, 
the mean values of the metallicity gradients are very similar 
in BCGs and ordinary ellipticals. 




-0.6 -0.4 -0.2 0.0 0.2 0.4 
log (r) arcsec 

Figure 2. Metallicity measurements normalised to an average 
mctallicity. 



In addition of doing 24 individual fits and averaging 
all these gradients, we also normalised the data to the same 
average metallicity value on the vertical scale (0.3) and then 
made a single fit to all the points simultaneously (see Figure 
[2]). Wc then find a mean gradient of -1.251±0.076 (with P 
= < 0.0001). 



Metallicity gradient - velocity dispersion relation 

A correlation between the metallicity gradient and the mass 
of galaxies is predicted by the classical dissipative col- 
lapse model (Larson 1974). Because dissipation efficiency 
increases with mass, we would expect to find stronger gra- 
dients in more massive galaxies. However, a relation be- 
tween the metallicity gradient and the mass of the galaxy 
- usually parametrised by the central velocity dispersion - 
has not been found in studies of normal elliptical galaxies 
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1.5 



1.0 - 




-1.0 - 



-1.5 J , , , , 1 

2.3 2.4 2.5 2.6 

log (a) km s" 1 

Figure 3. Metallicity gradients plotted against central velocity dispersion. The BCGs in our sample are shown with black symbols (those 
observed within 15 degrees of the MA (where known) are indicated with stars, and the others with circles. The BCGs from other studies 
are shown with empty symbols (those from Mehlert et al. (2003) with circles and those from Brough et al. (2007) with stars). The fitted 
metallicity - mass correlation for our sample is shown with a solid black line, and the correlation fitted to all the data with a dashed 
line. 



(for example Kuntschner et al. 2006; Sanchez-Blazquez et 
al. 2006b, 2007, although this relation may exists in clus- 
ter galaxies, see, for example, Forbes, Sanchez-Blazquez & 
Proctor 2005). A stronger metallicity gradient is found in 
intermediate luminosity galaxies while the relation between 
the metallicity gradient and the velocity dispersion gets flat- 
ter for more massive galaxies. The change in slope occurs 
at the transition between those galaxies supported by rota- 
tion and those galaxies dominated by random motions (see 
Sanchez-Blazquez et al. 2007). 

In this section we explore if BCGs behave differently 
from ordinary ellipticals. To do this, we plot the 24 metallic- 
ity gradients against central velocity dispersion (measured in 
the central a c /8 apertures) in Figure[3] We find a correlation 
between the two parameters with a 95 per cent confidence 
level (a t-test gives a probability of seven per cent that the 
two parameters are not correlated). A linear fit weighted by 
the errors in both, the velocity dispersion and the metallic- 
ity gradients gives a slope of -1.616T0.539. It can be seen 
that, contrary to the previous findings for elliptical galaxies 
- in the same mass range - not in the centres of clusters, the 
gradients seem to become steeper as the velocity dispersion 
increases. As an additional test, we also calculate the non- 
parametric Spearman rank coefficient (which is independent 
of the errors) as R s = -0.41 (for N =24). 

We also plot, in Figure[3] two BCGs from Mehlert et al. 
(2003; observed along the MA) and five from Brough et al. 
(2007; their sixth galaxy was excluded because it falls out 
of our mass range with a velocity dispersion of logo" = 2.18 
km s — l , and these galaxies were also observed along the 
major or intermediate axes) on the same graph. When we 
include these seven BCGs in the linear fit, regardless of slit 
orientation, we find a correlation consistent with a slope of 
-1.354 ± 0.427. Similarly, we find a probability of four per 



cent (P=0.0036) that these two parameters are not related 
according to a statistical t-test at 95 per cent confidence 
level, i.e. the correlation gets even stronger. 

Brought et al. (2007) also found, for their smaller sam- 
ple, a relation between the metallicity gradient and the 
velocity dispersion similar to the one observed here, but 
claimed it to be an artifact of the small size of the sample. 
However, we show here that the correlation persists even 
with a much larger sample of galaxies. 

As we have mentioned before, for normal elliptical 
galaxies, Davidge (1992) and Sanchez-Blazquez et al. (2007) 
found correlations between metallicity gradient and the 
anisotropj|^] We looked for this correlation in our sample 
of BCGs galaxies, using the V max and <to derived in Paper 1. 
Contrary to the works cited above, we do not find a correla- 
tion between these two parameters. Similarly to our study, 
Gorgas et al. (1990) also did not find such a correlation in a 
sample that combined BCGs and SOs. Whether these differ- 
ences are due to the presence or absence of BCGs in samples 
remain to be seen with larger samples. 

Metallicity gradients vs luminosity 

While the velocity dispersion is a proxy for the dynamical 
mass of the galaxy, the luminosity in the isT-band is sensitive 
to the predominantly red population in massive elliptical 

6 The amount of flattening due to rotation in a galaxy depends 
on the balance between ordered and random motions, and this 
can be quantified using the anisotropy parameter, which is de- 
fined as (Vmax/cro)* = ( V max /(ro )/ y/e/1 - t (Kormendy 1982). 
The rotational velocity (Vinax) is half the difference between the 
minimum and maximum peaks of the rotation curve, and ctq is 
the measured central velocity dispersion (as derived in Paper 1). 
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Figure 4. BCG metallicity gradients plotted against ii"-band lu- 
minosity. 



galaxies and hence, is very well correlated with the stellar 
mass (Brinchmann & Ellis 2000) . To check if the correlation 
found in the previous section is robust, we also check the re- 
lation between the metallicity gradients and the luminosity 
in the 7^-band. 

We plot the 2MASS absolute K-band luminosity 
against the metallicity gradients in Figure [4] We do a linear 
fit inversely weighted by the errors on the gradients for all 
the BCGs, regardless of slit orientation, and find a correla- 
tion with a slope of 0.f73 ± 0.081. We find a probability 
of five per cent (P=0.05) that these two parameters are 
not related according to a statistical t-test at 95 per cent 
confidence level. The non-parametric Spearman rank coeffi- 
cient (which is independent of the errors) is Rs = 0.31 (for 
N =22). We do not correct the .ff -magnitudes for passive 
evolution. Using the Bruzual & Chariot (2003) stellar popu- 
lation synthesis code with the assumption that the galaxies 
are 10 Gyr old and formed in an instantaneous burst, this 
correction is only -0.2 mag for a galaxy at z ~ 0.054 (in the 
if-band). Thus, it will not make a significant difference to 
whether or not a correlation is found. At the low redshift of 
this sample (mean z — 0.037 ± 0.003), the isT-corrections to 
the magnitudes are negligible. 



5.2 Age gradients 

From Table [3] it can be seen that we find only three signif- 
icant age gradients, and a mean gradient of 0.069 ± 0.049 
(and a standard deviation of 0.241), which is consistent with 
a zero age gradient as also found by Brough et al. (2007) for 
BCGs (0.01 ± 0.04). Mehlert et al. (2003) and Sanchez- 
Blazquez et al. (2007) also find null or very shallow age 
gradients for their samples of elliptical galaxies, and Spo- 
laor et al. (2010) also finds a mean age gradient of zero for 
their literature compilation of high mass early-type galax- 
ies. Fisher et al. (1995b) and Carter et al. (1998) found that 
their smaller samples of BCGs were younger in the centres 
than in their outer parts, but Cardiel et al. (1998a) found 
this only in cooling flow clusters. 



5.3 a-enhancement gradients 

In agreement with previous studies of both BCGs and nor- 
mal elliptical galaxies, we find null a-enhancement gradients 
in our sample of galaxies (with a mean gradient of -0.008 ± 
0.032 and with a standard deviation of 0.157), as shown in 
Figure [f]£] For individual galaxies, we only find one with an 
statistically significant gradient. 

We have used the well-calibrated indices H/3, Mgb 
and <Fe> to determine SSP-equivalent parameters, as it 
is known that stellar population parameters, in particular 
metallicity, depend strongly on the indices used to deter- 
mine it (Sanchez-Blazquez et al. 2006a, Pipino et al. 2010). 
However, we have tested whether using other indices affected 
the derived gradients. For example, we also used the indices 
H7A, H7F, C24668, Mgb with the same models to derive log 
age, [Z/H], and [E/Fe]. We then find zero significant age gra- 
dients, nine significant metallicity gradients, and one signif- 
icant a-enhancement gradient (mostly for the same galaxies 
for which we found significant gradients using the H/3, Mgb 
and <Fe> indices). Using H7A, H7F, C24668 and Mgb, we 
find a mean [Z/H] gradient of -0.289 ± 0.089, which is in 
agreement with -0.31 ± 0.05 found by Brough et al. (2007). 
It does not influence any of the results drawn from the SSP 
gradients in this study, as the gradients agree within the er- 
rors. We also did a boot-strap Monte Carlo simulation (with 
replacement) to assess the effect of non Gaussian errors on 
the significance quoted for the correlations. We found new 
correlations for [Z/H] gradients vs. velocity dispersion -2.184 
± 0.202 (instead of -1.616 ± 0.539) and [Z/H] gradients vs. 
K-magnitude 0.109 ± 0.015 (instead of 0.173 ± 0.081). 



5.4 Velocity dispersion gradients 

We also characterise the velocity dispersion gradients as 
A log(er)/A log(r) to enable a direct comparison with pre- 
vious studies, for example Brough et al. (2007). We fit the 
gradients using the same procedure as used for the indices 
and SSP parameters, and once again exclude the central 0.5 
arcsec (to each side) from the fit and take the gradient to be 
significant if it is larger than three times the la error on that 
gradient. The velocity dispersion gradients are tabulated in 
Table H 

We find negative velocity dispersion gradients for 17 
BCGs in this sample (of which six are bigger than three 
times the la error), and seven positive velocity dispersion 
gradients (of which two are bigger than three times the la 
error). In a much smaller sample of BCGs, Brough et al. 
(2007) found five out of their six velocity dispersion gradi- 
ents were consistent with being zero, and one negative gra- 
dient. As discussed in Paper 1, positive velocity dispersion 
gradients may imply rising mass-to-light (M/L) ratios with 
distance from the centre of the galaxies. 



7 The distribution of the a-enhancement gradients in Figure [5] 
can be non-symmetrical, and if the non-symmetrical lcr Gaus- 
sian error is read from the histogram then the error on the mean 
changes negligibly, and the mean is still zero within the errors 
(whether measured symmetrically or non-symmetrically). 
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Figure 5. Histograms indicating the mean log (age), a-enhanccmcnt, mctallicity and velocity dispersion gradients of our BCG sample. 



6 CLUSTER PROPERTIES 

The correlation between the metallicity gradient and the 
central velocity dispersion, which is absent in normal ellip- 
tical galaxies, suggests that the position of the galaxy in the 
centre of the cluster potential well may have a stronger influ- 
ence in the spatial distribution of line-strength indices than 
the mass of the galaxy. We explore the dependence of the 
stellar population gradients on the properties of the cluster, 
by investigating the following host cluster properties: X-ray 
luminosity Lx; X-ray temperature Tx; Cluster velocity dis- 
persion a duster', the BCG offset from the X-ray peak R s; 
and whether or not the host cluster is a cooling flow clus- 
ter or not. The literature values used here are also given in 
Table 7 in Paper 2, but are repeated here for completeness 
(Table |Eb. All the values are from spectra observed in the 
0.1 - 2.4 keV band with ROSAT, and using the same cos- 
mology, namely the Einstein de-Sitter model of Ho = 50 
km s" 1 Mpc" 1 , Qm — 1 and Qa — 0. As a test, the X-ray 
luminosity values were converted to the concordance cosmo- 
logical model (fi m = 0.3 and Q,a = 0.7) by calculating the 
appropriate cosmological luminosity distances of the clus- 
ters. We found that the assumed cosmological model does 
not influence the relative correlations between the parame- 
ters, and the previously most often used Einstein de-Sitter 
X-ray luminosity values are used to plot the correlations (as 
originally published from the ROSAT data). The other X- 



ray properties used here, such as X-ray temperature, do not 
depend on the cosmological model assumed (White, Jones 
& Forman 1997; Bohringer et al. 2004). The X-ray offset 
R g was calculated using Ho = 75 km s -1 Mpc^ 1 in Mpc. 
For the majority of the clusters this offset was not already 
published, and it was then calculated from the BCG and 
published X-ray peak coordinates. However, this was not 
possible for those clusters, e.g. Coma, where a BCG is not 
in the centre and where the coordinates of a corresponding 
local X-ray maximum were not available. 

Numerical simulations predict that the offset of the 
BCG from the peak of the cluster X-ray emission (R s) 
is an indication of how close the cluster is to the dynami- 
cal equilibrium state, and that this decreases as the cluster 
evolves (Katayama et al. 2003). We plot the metallicity gra- 
dients against R s in Figure [6] and do a linear fit inversely 
weighted by the errors for all the BCGs, regardless of slit ori- 
entation, to find a correlation with a slope of -7.546 ± 2.752. 
We find a probability of less than one per cent (P=0.009) 
that these two parameters are not related according to a sta- 
tistical t-test at 95 per cent confidence level. This agrees with 
our findings for the Mg2 gradients in Paper 4, and we find 
no other correlations with host cluster properties. The non- 
parametric Spearman rank coefficient (which is independent 
of the errors) is rather weak at Rs = —0.016 (for N =19). 
We again did a boot-strap Monte Carlo simulation (with 
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Table 4. BCG stellar velocity dispersion gradients. A * indicates 
where the gradient is bigger than three times the error on that 
gradient. 



replacement) to assess the effect of non Gaussian errors on 
the significance quoted for the correlations, and found the 
correlation: [Z/H] gradients vs. R ft -5.124 ± 0.696 (instead 
of -7.546 ± 2.752). 



7 DISCUSSION AND SUMMARY 

The luminosity and photometric uniformity of BCGs imply 
that they are not just the bright extension of the luminos- 
ity function of cluster galaxies (Bernstein & Bhavsar 2001). 
The position of these galaxies at the centres of clusters and 
their unique properties link their formation and evolution 
to that of their environment. However, the mechanisms be- 
hind their growth are still poorly understood. If environment 
has an effect on the star formation, the number of interac- 
tions a galaxy undergoes, or the dissipation of gas within the 
galaxy's potential well, we would expect to see an environ- 
mental dependence on the inferred stellar population gradi- 
ents. We have derived the gradients of the SSP-equivalent 
parameters (age, [Z/H] and [E/Fc]) and the velocity disper- 
sion for 24 BCGs in the nearby Universe. The host clusters 
of these galaxies comprise a large range of richness, ranging 
from groups to massive clusters. This sample is the largest 
BCG sample for which spatially resolved spectroscopy has 
been obtained. 

We find very shallow gradients in all three consid- 
ered parameters: age, metallicity and cv-enhancement. For 
normal elliptical galaxies, shallow gradients in age and a- 
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Figure 6. The upper plot shows the metallicity gradients against 
projected distance between BCG and X-ray peak of the host clus- 
ter. The lower plot shows the metallicity gradients against the 
stellar velocity dispersion of the galaxies where different symbols 
indicate whether or not the BCG is hosted by a cooling flow clus- 
ter (the grey stars indicate the cooling flow clusters, the empty 
circles indicate the non-cooling flow clusters, and the black filled 
circles indicate the clusters where the cooling flow status is not 
known). 

enhancement have also been found, and the mean metal- 
licity gradient is normally around —0.3 dex per decade of 
variation in the radius, similar to what we found here. We 
note that Gorgas et al. (1990) found that the Mg2 gradients 
of three BCGs are shallower than the mean gradient of nor- 
mal ellipticals. However, Davidge & Grinder (1995) found 
the D4000 gradients of six BCGs to be steeper than that 
for non-BCGs. The contradictory results could be the con- 
sequence of the limited sample sizes. We do not know if the 
large fraction of null metallicity gradients that we found is 
a consequence of the privileged position of the BCGs in the 
centre of the clusters or if it is due to their large masses. 

Up to now, there have been contradicting results among 
different studies about the existence of a relation between 
the metallicity gradients and the mass (usually param- 
eterised by the velocity dispersion) of elliptical galaxies 
(Peletier et al. 1990; Gorgas et al. 1990; Gonzalez 1993; 
Davies, Sadler & Peletier 1993; Carollo & Danziger 1994; 
Mehlert et al. 2003; Forbes et al. 2005), although a con- 
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Table 5. X-ray properties and velocity dispersions of the host clusters. The cluster velocity dispersion (<T c i ustor ) values are in km s _1 
and the projected distance between the galaxy and the cluster X-ray peak (i? ff ) is in Mpc. The presence of a cooling flow is indicated 
with a /, the absence of a cooling flow with a X, and - indicates that the cooling flow status is not known. The * marks at R ff indicate 
that the galaxy is not in the centre of the cluster but closer to a local maximum X-ray density, different from the X-ray coordinates 
given in the literature. The references are: a = Chen et al. (2007); b = Bohringer et al. (2004); c = Bohringcr ct al. (2000); d = White 
(2000); f = Peterson ct al. (2003); g = Sivanandam ct al. (2009); h = David ct al. (2003); i = Vikhlinin et al. (2009); j = Osmond & 
Ponman (2004); k = Sakelliou, Hardcastle & Jetha (2008); 1 = Finoguenov, Arnaud & David (2001); n = Ikebe et al. (2002); o = Sun ct 
al. (2009); q = Cavagnolo ct al. (2008); r = Fukazawa, Makishima & Ohashi (2004); w = White, Jones & Forman (1997); c = Edwards 
ct al. (2007); cc = Calculated from Bohringcr ct al. (2000); cl = Calculated from Ledlow et al. (2003); cb = Calculated from Bohringcr 
ct al. (2004); m = Mahdavi & Gcllcr (2001); st = Struble & Rood (1999); s = Struble & Rood (1991); p = Patel et al. (2006); pe = 
Peres et al. (1998). All X-ray luminosity (Lx) values were measured in the 0.1-2.4 keV band. 



sensus started to appear that the the behaviour is differ- 
ent in different mass ranges (see Carollo, Danziger & Buson 
1993; Sanchez-Blazquez et al. 2007; Spolaor et al. 2009; and 
Kuntschner et al. 2010). For elliptical galaxies, gradients get 
steeper with mass for galaxies with masses below ~10 n M Q , 
while the opposite behaviour is observed for more massive 
galaxies (although the slope is much flatter and the scatter 
also increases). We find a correlation with a slope of -1.616 
± 0.539 (P=0.007) for the metallicity gradients against cen- 
tral velocity dispersion of the 24 galaxies in this sample. 
This correlation becomes slightly stronger when we also in- 
clude the two BCGs from Mehlert et al. (2003) and five from 
Brough et al. (2007) in the same mass range. Contrary to 
the previous findings for elliptical galaxies not in the cen- 
tres of clusters, the gradients seem to become steeper as the 
velocity dispersion increases. 

We notice that this steepening of BCG gradients was 
also shown in Spolaor et al. (2009) using a literature com- 
pilation of eight BCGs. We also find a correlation between 
the metallicity gradients and K-b&nd luminosity supporting 
that there is a trend between the strength of the gradient 
and the mass of the BCGs that it is not present, or at least 
not so evident in elliptical galaxies. 



To study the possible influence of the cluster on the 
properties of the central galaxy, we investigated the metal- 
licity gradients against the host cluster properties (X-ray 
luminosity Lx; X-ray temperature Tx; cluster velocity dis- 
persion a c i us ter', the BCG offset from the X-ray peak R g; 
and whether or not the host cluster is a cooling flow cluster 
or not). We find a strong correlation between the metallicity 
gradients and R ofi with a slope of -7.546 ± 2.752 (P=0.009). 
This agrees with our findings for the Mg2 gradients in Paper 
4. We find no other correlations with host cluster properties. 
Notice that Cardiel et al. (1998a) found differences in the 
Mg2 gradients between galaxies in clusters with and without 
cooling flows. 



In summary, we have found hints of differences in the 
spatially resolved stellar population properties of BCGs and 
normal elliptical galaxies pointing to an influence of the clus- 
ter in the way these systems assemble. The trends showed 
here are, admittedly, weak, but they are interesting as this 
is the first time they have been found. These correlations 
are worth exploring in more detail, with stellar population 
gradients reaching larger radii. 
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APPENDIX A: COMPARISON WITH 
PREVIOUS MEASUREMENTS 

To test self-consistency, we compared the indices measured 
in central bins with the central measurements (from Paper 



2, measured in a c /8 apertures) in Figure Al 

The index gradients of one of the galaxies could also be 
compared to previous measurements in the literature from 
Fisher et al. (1995b) and Mehlert et al. (2000), as shown in 
Figure |A2| The profiles compare well within the errors and 
small differences can be as a result of spectral resolution and 
slit width. 



APPENDIX B: INDEX AND STELLAR 
POPULATION GRADIENTS 



The H/3 index gradients are plotted in Figures |B1| and |B2| 
and the Mg6, Fe5270 and Fe5335 index gradients are plotted 
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Figure Al. Comparison between the central index measurements from Paper 2 (measured in a e /8 apertures) and the central bins 
measured in this work (minimum S/N 40 per A at H/3) for four key indices. 
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Figure A2. Comparison to previous data (NGC4839). The data are folded with respect to the centre of the galaxy. Mg2 is measured in 
magnitude, and all the other indices shown are measured in A. 
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Figure Bl. BCG H/3 gradients. All data is folded with respect to the galaxy centres. The central 0.5 arcsec, to each side (thus in total 
the central 1 arcsec, comparable to the seeing), were excluded in the figure as well as in the fitted correlations. 
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Figure B2. BCG H/3 gradients continue. 
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Figure B3. BCG Mg;,, Fe5270 and Fe5335 gradients. All data is folded with respect to the galaxy centres. The central 0.5 arcsec, to 
each side (thus in total the central 1 arcsec, comparable to the seeing), were excluded in the figure as well as in the fitted correlations. 
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Figure B4. BCG Mg 6 , Fc5270 and Fe5335 gradients continue. 
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Figure B5. BCG Mg;,, Fc5270 and Fe5335 gradients continue. 
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Figure B6. BCG Mg b , Fe5270 and Fe5335 gradients continue. 
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Figure B7. Age, a-enhanccment, metallicity and velocity dispersion profiles of the BCGs. The galaxies were binned in the spatial 
direction to a minimum S/N of 40 per A in the H/3 region of the spectrum. The central 0.5 arcsec, to each side, are excluded. In all the 
gradient plots, the data were folded with respect to the galaxy centres. 
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Figure B8. Age, a-enhancement, metallicity and velocity dispersion profiles of the BCGs continue. 
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Figure B9. Age, ct-enhancement, metallicity and velocity dispersion profiles of the BCGs continue. 
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Figure BIO. Age, a-enhancement, metallicity and velocity dispersion profiles of the BCGs continue. 



